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■ Abstract 

^■f^ ' Further analysis has been made on the miUimcter and submilHmeter-wave (100-1600 GHz or 3 mm - 

I 188 /im) atmospheric opacity data taken with the Fourier Transform Spectrometer (FTS) at Pampa la Bola, 

• 4800 m above sea level in northern Chile, which is the site of the Atacama Large Millimetcr/submillimeter 

^ I Array (ALMA) . Time-sequence plots of millimeter and submillimeter-wave opacities show similar variations 

to each other, except for during the periods with liquid water (fog or clouds) in the atmosphere. Using 
millimeter and submillimeter-wave opacity correlations under two conditions, which are affected and not 
affected by liquid water, we succeeded to separate the measured opacity into water vapor and liquid water 
' opacity components. The water vapor opacity shows good correlation with the 183 GHz water vapor 

, line opacity, which is also covered in the measured spectra. On the other hand, the liquid water opacity 

■ and the 183 GHz line opacity show no correlation. Since only the water vapor component is expected to 
I affect the phase of interferometers significantly, and the submillimeter-wave opacity is less affected by the 

liquid water component, it may be possible to use the submillimeter-wave opacity for a phase-correction 
I ' of submillimeter interferometers. 

2 i Key words: atmospheric effects — instrumentation: spectrometer — site testing — submillimeter 

+13 . 

c/2 . 

■ 1. Introduction Chamberlin 2000: Radford et al. 2001; Radford 2002; 
J> , Sakamoto 2002), a 350 pmi tipping photometer (Radford, 

In astronomy, observations with millimeter and sub- Holdaway 1998; Radford 2002), 183 GHz water line ra- 

^ , millimeter wavelengths are highly suitable for study- diometers (Delgado et al. 1999, 2000, 2001), and 11 GHz 

■ ing molecular gas and dust emission from various radio seeing monitors (Radford et al. 1996; Holdaway et 
astronomical objects, including young stellar objects, al. 1997; Radford, Holdaway 1998; Butler et al. 2001; 
star-forming/starburst regions, and primeval galaxies. Radford 2002; Sakamoto 2002) are now continuously op- 
However, observations at these frequencies/ wavelengths, erating to evaluate the absorption and phase fluctuation 
especially at submillimeter-wavelength, are difficult, since caused by water constituents. Also, a 492 GHz tipping ra- 
it is easily and heavily absorbed by water constituents diometer performed measurements for intermittent peri- 
in the atmosphere. Therefore, understanding the absorp- ods of time (Hirota et al. 1998). We then carried out wide 
tion/emission spectra caused by the water constituents is frequency coverage opacity measurements using a Fourier 
important, but is not yet well characterized. Transform Spectrometer (FTS) in the winter of 1997 and 

At the high-altitude site of the Atacama desert in 1998. 
Northern Chile (Pampa la Bola - Chajnantor region. Opacity measurements with the FTS were performed 

about 5000 m above the sea level), the US, Europe, during 1997 September 5 and 12. Opacity correlations be- 

and Japan are planning to build the Atacama Large tween 220 GHz and all the submillimeter-wave windows 

Millimeter/submillimeter Array (ALMA), which will be were obtained, although the weather conditions were lim- 

the world's largest millimeter and submillimeter interfer- ited (Matsuo et al. 1998a, b). During 1998 June 16 and 

omcter. Since the site is very dry, it is a suitable loca- 18, opacities under good weather conditions were mea- 

tion for studying the behavior of water constituents in sured, and high transmission spectra as well as the detec- 

the millimeter/submillimeter region. At Pampa la Bola - tion of supra-terahertz windows located around 1035 GHz, 

Chajnantor area, 220 GHz/225 GHz tipping radiometers 1350 GHz, and 1500 GHz have been obtained (Matsushita 

(Kohno et al. 1995; Radford, Holdaway 1998; Radford, et al. 1999, 2000; Matsushita, Matsuo 2002). These 
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supra-terahertz windows have also been detected by other 
groups (Paine et al. 2000; Pardo et al. 2001b, 2002) and 
successfully modeled with radiative transfer calculations 
(Pardo et al. 2001a,b, 2002; Matsushita et al. 1999, 2000). 
Combining these data sets with long-term site testing ra- 
diometer data (Radford, Holdaway 1998), the fraction 
of time with a low-opacity condition in submillimeter- 
wave (675 GHz and 875 GHz) atmospheric windows at 
the ALMA site was estimated. The submillimeter zenith 
opacity is less than 1.0 for about 50% of the winter season 
and about 37% of a year, which is about 50% better than 
the Mauna Kea site (Matsushita et al. 1999). 

During the measurement period in 1998, quite different 
opacity correlations compared with the normal ones men- 
tioned above were obtained. Because of their behavior, 
weather condition, and theoretical perspectives, it is sug- 
gested that liquid water in the atmosphere caused these 
correlation changes (Matsushita et al. 1999; Matsushita, 
Matsuo 2002). We analyzed these 'abnormal' opacity cor- 
relation data in detail, and considered the effect of liquid 
water on millimeter- and submillimeter- waves. Based on 
these results, we propose new phase-correction methods 
for submillimeter interferometry. 

2. Measurements and Data Analysis 

Atmospheric opacity measurements were performed at 
Pampa la Bola (northern Ghile, the Atacama desert, 
4800 m altitude) using a Martin-Puplett type FTS 
(Martin, Puplett 1970) and an InSb bolometer. The fre- 
quency coverage was 100-1600 GHz (or 3 mm to 188 /im 
in the wavelength range). Two different observing runs 
were carried out, the first one in 1997 September 5-12, 
and the second one in 1998 June 16-18 (both in winter 
season in Chile). 

The beam size of the instrument was about 10°, and 
the apodized frequency resolution was ~ 0.3 cm^^ (^ 
10 GHz). Atmospheric emission spectra were obtained to- 
ward different air masses [sec(z), z is zenith angle] at 1.0, 
1.5, 2.0, 2.5, and 3.0 by rotating a tipping mirror outside of 
the spectrometer. A tipping measurement with the FTS 
including two temperature calibrations was made every 
12 min in the first observing run and every 14 min in the 
second run. The absolute-brightness temperature of the 
atmosphere was calibrated using blackbodies (Eccosorb 
AN74) at liquid-nitrogen temperature (~ 73 K at 4800 m 
altitude) and at ambient temperature. The ambient tem- 
perature was monitored with a weather station (Sakamoto 
et al. 2000). In data analysis, tipping measurements were 
used for frequencies lower than 450 GHz, and tempera- 
ture measurements were used for frequencies higher than 
450 GHz. From 10 PM on June 17 to 4 AM on June 18 
in 1998, however, temperature measurements were used 
for all frequencies. This is because the existence of clouds 
in the atmosphere is suggested (see subsection 4.1), and 
therefore requirements for tipping measurement (need a 
constant sky condition toward different air masses and a 
low opacity condition) would not be satisfied. Further 
details on the instruments, measurements, calibrations, 



and data analysis are described elsewhere (Matsuo et 
al. 1998a,b; Matsushita et al. 1999). 

Since the frequency coverage of our FTS data ranges 
from 100 GHz, we can measure the 183 GHz water va- 
por line opacity by subtracting the continuum component 
from both sides of the line. We assumed 161-204 GHz 
as the 183 GHz water vapor line frequency range, and 
136-161 GHz and 204-230 GHz as contimmm frequency 
ranges. The continuum opacity within the line frequency 
range was calculated with the interpolation of the contin- 
uum opacities located at both sides of the line. After con- 
tinuum subtraction, the line opacity was divided by the 
number of the data points within the line frequency range 
(12 points), since our FTS frequency resolution is not good 
enough to resolve the 183 GHz line. The 183 GHz opacity 
used in this paper is, therefore, the average opacity in this 
frequency range. 

3. Time- Variations of Measured Opacities 

Time- variation plots for 220, 345, 492, 675, 875, and 
937 GHz opacities are shown in figure 1. The selected 
frequencies are either located at the center of atmospheric 
windows and/or are frequencies of special astronomical 
interests (see Matsushita et al. 1999; Matsushita, Matsuo 
2002). The left-hand and right-hand sides of the figure are 
the data taken in 1997 and 1998, respectively. In the 1997 
measurement, because the sensitivity of the bolometer in 
daytime data had large variations, we only used nighttime 
data. 

The time- variation plots of the 1997 data (the left-hand 
side of figure 1) clearly show that the opacity condition 
was relatively good at the beginning of the measurement 
(1997 September 5), but gradually worsened toward the 
end (1997 September 12). The overall time variations of 
the plots look similar to each other, except for the plot of 
937 GHz which saturates under worse weather conditions. 

On the other hand, the time-variation plots of the 1998 
data (the right-hand side of figure 1) indicate that the 
opacity condition was extremely good during the first 
half of the measurement (from the evening of June 16 
to around noon of 17), but gradually worsened from the 
afternoon of June 17, and became totally bad around mid- 
night. The lowest opacity (the best atmospheric transmis- 
sion; see figure 1 in Matsushita et al. 1999) was recorded 
around noon on June 17 (Chilean local day of around 
17.5 in the plots). The first half of the time- variation 
plots looks similar for all frequencies (from 220 GHz to 
937 GHz), but the second half, especially around the be- 
ginning of June 18, of the plots is different between the 
lower frequencies (220 GHz and 345 GHz) and higher fre- 
quencies. The implications of this is discussed in a later 
section. 

Figure 2 shows time-variation plots of the supra- 
terahertz windows. Since the weather condition of the 
1997 measurement was not good enough to detect these 
supra-terahertz windows, we only show the 1998 data. 
Because the sensitivity of our FTS around the 1350 GHz 
and 1500 GHz windows is low and has some systematic 
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Fig. 1. Time-sequence plots for 220, 345, 492, 675, 875, and 937 GHz opacities. Tlie left-hand and right-hand sides of the figure 
arc the data taicen in 1997 September and 1998 June, respectively. The Chilean local time corresponds to UT — 4 hr. 
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Fig. 1. (Continued) 
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Fig. 2. Time-sequence plots for (a) 1035 GHz, (b) 1350 GHz, 
and (c) 1500 GHz opacities. The Chilean local time corre- 
sponds to UT — 4 hr. Since the weather conditions during 
the 1997 measurement were not good enough to detect these 
supra-torahertz windows, we only show 1998 data plots. 
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Fig. 3. Time-sequence plots for 183 GHz water vapor line 
opacity averaged over 161 to 204 GHz, taken in (a) 1997 
September and (b) 1998 June. The Chilean local time corre- 
sponds to UT - 4 hr. 

errors (see Matsushita ct al. 1999 for detail), the plots 
have a larger scatter than the other plots. Still, the plots 
show a similar trend to the others. 

We also plot the time variation of the 183 GHz water va- 
por line opacity (figure 3). The 1997 data plot (figure 3a) 
shows similarities to the other opacity plots (the left-hand 
side of figure 1). The 1998 data plot (figure 3b) shows sim- 
ilar variation to the higher frequency (> 492 GHz) plots 
(the right-hand side of figure 1 and figure 2), but different 
from the lower frequency (220 GHz and 345 GHz) plots, 
especially from the late evening of June 17 to the early 
morning of June 18. 

4. Analysis of Water Vapor and Liquid Water 
Opacity Components 

4-.1. 'Abnormal' Opacity Correlations between 
Millimeter- and Submillimeter- Wave 

Usually, correlation plots between 220 GHz opacity and 
submillimeter window opacities are strongly correlated 
with coefficients of 22-24 (hereafter we refer to these cor- 



6 



S. Matsushita and H. Matsuo 



relations as 'normal' correlations; Matsuo et al. 1998a; 
Matsushita et al. 1999, 2000; see also Masson 1994). 
However, from 17.9 to 18.2 Chilean local day (from 10 PM 
on June 17 to 4 AM on June 18) in 1998, we obtained quite 
different correlation plots from the 'normal' correlations 
(see figure 4 in Matsushita et al. 1999). These 'abnor- 
mal' plots show lower correlation coefficients and larger 
zero-point offsets than those of the 'normal' correlations, 
which means that the fractional variation of the millimeter 
opacity is much larger than the submillimeter opacity. 

The NRAO surveillance camera at the Chajnantor site 
(http://www.tuc.nrao.edu/alma/site/site.html) shows a 
clear sky image around 5 PM on June 17, but shows a 
cloudy image around 9 AM on June 18. In addition, the 
ambient temperature data taken by the weather station 
(Sakamoto et al. 2000) recorded temperature below freez- 
ing (around — to — 8°C). These weather data suggest 
that there were water droplets (liquid water), or ice parti- 
cles, in the atmosphere when we obtained the 'abnormal' 
correlation data. 

Water vapor absorbs stronger in the submillimeter-wave 
region than in the millimeter-wave region (Liebe 1989; 
Liebe et al. 1993; Pardo et al. 2001a). Liquid water (cloud 
or fog); on the other hand, absorbs more effectively at 
the millimetcr-wavc region for a given particle-size dis- 
tribution (Ray 1972; Liebe 1989; Liebe et al. 1991; see 
also Matsushita et al. 2002), which means that a vari- 
ation in the quantity of the liquid water content in the 
atmosphere causes a large change in the millimcter-wave 
opacity. The existence of ice cannot explain the 'abnor- 
mal' correlations. Ice has about two orders of magni- 
tude lower values for the imaginary part of the index of 
refraction (« absorption coefficient^) at millimeter- wave 
region than that of liquid water (Warren 1984; Liebe et 
al. 1993; Hufford 1991; Ray 1972), and is similar to that 
of water vapor. 

There is a possibility that other molecules or aerosols 
cause the 'abnormal' correlations. Neither the oxygen 
nor ozone lines can be the cause, since these show up 
with narrow lines in the spectra (e.g., Pardo et al. 2001a; 
Matsushita et al. 1999). The characteristics of the wet and 
dry pseudocontinuum components, which are thought to 
be the far wings of water vapor lines in the infrared range 
and collision-induced absorption of the dry atmosphere 
(N2 and O2), respectively (Pardo et al. 2001a,b, 2002), are 
still empirical and not well understood. Furthermore, the 
wet pseudocontinuum should be correlated with water va- 
por, and the dry atmosphere should have a constant den- 
sity and is not expected to fluctuate. In addition, since the 
wet and dry pseudocontinuum absorptions at millimeter- 
wave are too weak to explain the 'abnormal' correlations, 
we have to introduce a new kind of pseudocontinuum ab- 
sorption component, which has rather stronger absorption 
in the millimeter- wave range than the other pseudocontin- 
uum components. However, it is beyond the scope of our 
paper to introduce the new pseudocontinuum component. 

^ The imaginary part of the complex index of refraction, nj, 
and the absorption coefficient, a, is related as a = Anuni/c 
(Thompson et al. 2001). 



Table 1. Correlation coefficients of the 'normal' 
(awv) s-nd 'abnormal' (clq) correlations between 
220 GHz opacity and submillimeter opacities. 



Frequency 


awv* 




(GHz) 






345 


3.65 ±0.02 


1.42 ±0.03 


410 


7.53 ±0.05 


1.64 ±0.05 


492 


23.6 ±0.3 


1.5±0.2 


675 


22.4±0.2 


1.9±0.2 


875 


24.2 ±0.2 


1.8±0.2 


937 


43.9 ±0.9 


0.6 ±0.2 



* The values are taken from Matsushita et al. (2000). 

We therefore suggest that the 'abnormal' correlations 
are caused by the effect of the liquid water absorption, 
which is the most plausible known component. The reason 
why liquid water still existed in the atmosphere even when 
the ambient temperature was below freezing may be that 
the atmospheric temperature at the clouds was warmer 
than that on the ground, or that the clouds consisted of 
both liquid water and ice, but the absorption caused by 
liquid water dominated the total atmospheric absorption. 

4-2. Separation of Total Opacity into Two Opacity 
Components 

To confirm this suggestion, we separate the total opacity 
into water vapor opacity and liquid water opacity compo- 
nents using the 'normal' and 'abnormal' correlations be- 
tween millimeter and submillimeter opacities. For the sep- 
aration, we assume that (a) the 'normal' correlations re- 
sult only from the water vapor absorption, (b) 'abnormal' 
correlations are the combined effect of the water vapor and 
liquid water absorptions, (c) the correlation coefficients of 
the liquid water absorption correspond to those of the 
'abnormal' correlations, and (d) the zero-point offsets are 
zero in all the correlations. Under these assumptions, we 



can write: 

Tu,WV = CLWV ■ T220,WV, (1) 

VLQ = ^LQ ■ '''220,LQ, (2) 

'''220 ~T 220, WV + T'220,LQ , (3) 

T,y = r^,wv + Ti,^LQ, (4) 



where T220 and t^, indicate opacities at 220 GHz and fre- 
quency v, and the suffix WV and LQ indicate water vapor 
and liquid water components, respectively. The known 
(measured) parameters are awvj Qlq, ''"220, and r^. Using 
these parameters, we can separate the water vapor and 
liquid water opacity components from the measured to- 
tal opacity. We summarized the awv and clq of the se- 
lected frequencies derived from the 'normal' (Matsushita 
et al. 2000) and 'abnormal' (Matsushita et al. 1999) cor- 
relations in table 1. 

We then performed the component separation on the 
220 GHz and 675 GHz opacities taken in 1998, and the re- 
sults of the separated millimeter and submillimeter opac- 
ities are shown in figure 4. The top row of figure 4 shows 
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Fig. 4. Time variations of (a) total 220 GHz opacity, (b) 220 GHz separated water vapor opacity (filled circles, see left-hand axis 
for their unit), overplotted with the same plot but multiplied by four (open circles, see right-hand axis), (c) 220 GHz separated 
liquid water opacity, (d) total 675 GHz opacity, (e) 675 GHz separated water vapor opacity, and (f) 675 GHz separated liquid water 
opacity. The data were taken during the 1998 measurement, and the Chilean local time corresponds to UT — 4 hr. 



the time variations at 220 GHz for (a) the total opac- 
ity, (b) the separated water vapor opacity, and (c) the 
separated liquid water opacity. The bottom row shows 
the same plots but for 675 GHz. The total opacity plots 
for 220 GHz and 675 GHz (figures 4a,d) look different, 
especially around Chilean local day from 17.9 to 18.2, 
but the separated water vapor opacity plots (figures 4b, e) 
show similar time variation to each other. On the other 
hand, the 220 GHz separated liquid water opacity plot 
(figure 4c) shows a quite different behavior from that of 
675 GHz (figure 4f), when the 'abnormal' correlation oc- 
curred. These results suggest that millimeter opacity is 
largely affected by the liquid water component in the at- 
mosphere, although the effect in submillimeter opacity is 
small. 

^.3. Correlations with 183 GHz Water Vapor Line 
Opacity 

To further confirm that the separated water vapor and 
liquid water opacities were really caused by these compo- 
nents in the atmosphere, it is very important to compare 
these with the 183 GHz water vapor line opacity. To this 
end, we made correlation diagrams between the 183 GHz 
line opacity and those of total, water vapor, and liquid wa- 
ter at 220 GHz (millimeter) and 675 GHz (submillimeter), 
as shown in figure 5. The correlation diagram between 
the 183 GHz line opacity and 220 GHz total opacity (fig- 
ure 5a) shows a large scatter, but after the separation, the 
220 GHz water vapor opacity shows good correlation with 
the 183 GHz line opacity (figure 5b). On the other hand, 
the 220 GHz liquid water opacity is completely indepen- 



dent from the 183 GHz line opacity (figure 5c). These 
results confirm the success of the opacity component sep- 
aration. 

The correlation diagram of the 183 GHz line opacity 
and the 675 GHz total opacity (figure 5d) is similar to 
that of the 183 GHz line opacity and the 675 GHz water 
vapor opacity. The correlation diagram of the 183 GHz 
line opacity and the 675 GHz liquid water opacity, on the 
other hand, shows little contribution of liquid water in the 
675 GHz opacity. These results support that the submil- 
limeter opacity is dominated by water vapor absorption, 
and less affected by liquid water absorption. 

The correlation diagrams between the water vapor 
opacity and the 183 GHz line opacity (figures 5b, e) show 
a somewhat curved structure. The curve at the lower end 
of the opacity, which shows an offset from the zero opacity 
of the water vapor opacity, is caused by the limitation of 
the measurements (see Matsuo et al. 1998a, b) and/or by 
both the wet and dry pseudocontinuum absorptions. The 
curve at the upper end can be explained by the satura- 
tion of the 183 GHz line. This saturation is theoretically 
expected (Lay 1998; Carilli et al. 1998) and also observed 
with multi-channel 183 GHz line radiometers (Wiedner, 
Hills 2000; Yun, Wiedner 1999; Wiedner et al. 2001). 

5. Application for Phase-Correction 

The most unwanted atmospheric effect for millime- 
ter/submillimeter interferometers is excess paths in elec- 
tromagnetic wave propagations caused by water vapor. A 
fluctuation in the path difference between a pair of anten- 
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Fig. 5. Correlation diagrams between 183 GHz water-vapor line opacity averaged over 161 to 204 GHz and (a) total 220 GHz 
opacity, (b) 220 GHz separated water vapor opacity, (c) 220 GHz separated liquid water opacity, (d) total 675 GHz opacity, (e) 
675 GHz separated water vapor opacity, and (f) 675 GHz separated liquid water opacity. The filled and open circles indicate the 
dataset taken in 1997 September and 1998 June, respectively. 



nas causes phase fluctuation, which decreases the inter- 
ferometric data quality noticeably. Since the characteris- 
tics of water vapor are theoretically and observationally 
well known, it is better to use its absorption for phase- 
correction methods. Liquid water, on the other hand, has 
a different index of refraction from that of water vapor 
(Thompson ct al. 2001; Ray 1972; Liebe 1989), and its con- 
tribution to the opacity at lower frequency is much higher 
than that of water vapor, but affects the phase fluctua- 
tions less (Matsushita et al. 2002). In addition, the char- 
acteristics of liquid water are poorly understood. Hence, 
for phase-corrections using opacity (brightness tempera- 
ture) of the atmosphere, one should measure the water 
vapor opacity instead of the opacity affected by liquid wa- 
ter. 

Based on the results and discussion in the previous sec- 
tion, wc suggest that an accurate phase-correction for in- 
terferometers can be made using (a) 183 GHz water vapor 
line radiometers, which avoids the line center that can 
be easily saturated, (b) submillimetcr-wave radiometers 
whose frequency range is much less affected by the liquid 
water absorption, or (c) two frequency (millimeter- and 
submillimeter-wave) radiometers which can separate the 
water vapor and liquid water components, as mentioned 
above. 

The phase-correction method (a) has already been dis- 
cussed and demonstrated by many authors (Wiedner et 
al. 2001, 2002; Delgado et al. 2001; Wiedner, Hills 2000). 

Method (b) needs either one radiometer in each an- 
tenna, or it is also possible to use the total power output 
from a submillimeter receiver for astronomical observa- 



tions, and hence no additional instruments are required. 
Although the effect of liquid water in the submillimeter- 
wave is much less than that in millimeter- wave, this 
method may be inaccurate under a large liquid water con- 
tent. We can estimate the phase-correction errors due 
to the liquid water content, assuming that it contributes 
to a typical 220 GHz opacity scatter of 10%. The esti- 
mated errors are 19 /im, 42 /im, and 71 /im under 220 
GHz opacities {T220) of 0.016 (the best condition at the 
ALMA site; Matsushita et al. 1999), 0.036, and 0.061 
(25% and 50% quartile of T220, respectively; Radford, 
Chamber lin 2000), respectively. These errors correspond 
to phase errors of 15°, 34°, and 57° at 675 GHz, respec- 
tively, which are the accuracies for this phase-correction 
method (see Matsushita et al. 2002 for detail). 

Method (c) needs two radiometers in each antenna or 
uses the total power outputs of both millimeter and sub- 
millimeter receivers that operate simultaneously. In spite 
of the system complexity, this method is advantageous, 
since there arc fewer errors due to the liquid water content. 
One important question is whether olq is constant or not. 
The emission properties of liquid water as a function of 
frequency will depend on droplet size and temperature, 
which are likely to change from night to night. More mea- 
surements under large liquid water content are needed. 
However, the general trend, i.e., submillimeter measure- 
ments are less affected by liquid water, does not change. If 
the uncertainty of the liquid water model is 10%, the error 
of the 675 GHz total power phase-correction mentioned 
above will be reduced by 1/10 (Matsushita et al. 2002). 

The merit of method (b) is that the field of view for 
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the phase-correction is exactly the same as that of the 
astronomical observations. A disadvantage of submillime- 
ter phase-correction methods is that they cannot be used 
to observe in submillimcter- wavelength under bad weather 
conditions, and therefore cannot be applied to a site where 
the submillimcter opacity is large. 

On the other hand, in some cases, the effect of ice parti- 
cles might exist. Since the absorption characteristics of ice 
particles are similar to those of water vapor (see subsec- 
tion 4.1), it may be difficult to distinguish between these 
two absorptions by the methods (b) and (c). In addition, 
it is not clear to which amount of ice particles usually exist 
in the atmosphere. These might be one of the uncertain- 
ties in the proposed total-power phase-correction meth- 
ods. To clarify these problems, further measurements of 
atmospheric opacity spectra due to ice particles and de- 
tailed theoretical modelings will be needed. 

To check the feasibility of the newly proposed meth- 
ods (b) and (c), test measurements with other phase- 
correction methods, such as the method (a), or the fast- 
switching method (Holdaway 1992; Holdaway, Owen 1995; 
CariUi et al. 1996; CariUi, Holdaway 1997, 1999; Morita et 
al. 2000) should be performed. The best available instru- 
ment for such tests is the Sub-Millimeter Array (SMA; 
Ho 2000) on Mauna Kca, Hawaii, which is being con- 
structed by the Smithsonian Astrophysical Observatory 
(SAO) and the Acadcmia Sinica Institute of Astronomy 
and Astrophysics (ASIAA) of Taiwan. SMA has a similar 
observing frequency range to ALMA (180-900 GHz), and 
will have a capability of dual-frequency operation. In ad- 
dition, two 183 GHz radiometers (Wiedner et al. 2001) 
are already installed, and improved 183 GHz radiome- 
ters (Hills et al. 2001) are being developed. SMA will 
therefore be an ideal instrument to evaluate the proposed 
total-power phase-correction methods and other various 
methods. 

6. Conclusion 

We have made further analyses on the FTS measure- 
ments of millimeter and submillimcter atmospheric opac- 
ity taken at Pampa la Bola. Part of the 1998 opacity 
data is heavily affected by liquid water in the atmosphere, 
which served to understand the absorption spectrum of 
liquid water. We successfully separated the measured 
opacity into water vapor and liquid water opacity compo- 
nents. The separated water vapor opacity shows a good 
correlation with the 183 GHz water vapor line opacity, 
and the liquid water opacity shows no correlation with 
the line opacity. 

The results also indicate that the submillimcter opacity 
is less affected by liquid water than the millimeter opac- 
ity; we therefore propose to use the submillimeter opac- 
ity for the interferometric phase-correction by using the 
total-power outputs from submillimeter receivers for as- 
tronomical observations. 
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